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摘 要 : 提出 -种 简单 而 有 效 的 参数 不 确定 时 灌 系 统 和 鲁 棒 PID 控制 器 设计 方法 . 通过 在 ,二 ; 平面 上 绘制 稳定 边界 
线 ,确定 稳定 的 PID 控制 器 参数 区 域 ;推导 了 一 阶 不 稳定 时 滞 系 统 PI 控制 器 和 PID 控制 器 的 存在 性 条 件 ;基于 推广 
到 时 灌 系 统 的 校 边 定理 ,确定 所 有 重 棒 PID 控制 器 参数 集 . 仿真 实例 表明 了 该 方法 的 优越 性 . 
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Robust control design for time-delay systems with parameter 
uncertainties using the PID controller 
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Abstract: A simple and efficient method for the design of robust PID controllers for time-delay systems with 
parameter uncertainties is presented. The method is based on plotting the stability boundary locus in the 大 plane 
and then calculating stabilizing values of the parameters of a PID controller. Moreover, a necessary condition for the 
exist ence of a PI controller or a PID controller to simultaneously stabilize the open-doop unstable time -delay systems 


is also derived. The complete set of robust PID controller parameters is determined by the edge theorem extended to 
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time-delay systems. Examples are given to show the benefit of the presented method. 
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引 写 

工业 控制 过 程 由 于 机 理 复杂 以 及 时 变 、 时 灌 、 
非 线性 和 耦合 等 原因 , 使 其 精确 的 数学 模型 很 难 建 
过 ,而 其 广义 传递 画 数 可 用 一 阶 惯 性 环节 加 一 纯 时 
灌 过 程 来 近似 ", 实践 证 明 这 种 近似 是 有 其 适应 性 
的 . 尽管 现代 最 优 控 制 技术 发 展 迅速 ,如 万 ;和 有 ~ 
鲁 棒 控制 技术 ”1 ,但 至 少 存在 以 下 两 方面 原因 , 使 
万: 和 万 ~ 控制 技术 并 不 能 直接 应 用 于 工程 实际 :一 
方面 由 于 扰动 的 引入 将 带 来 控制 器 设计 的 保守 性 ; 
男 一 方面 由 于 石和 矿 ~ 控制 器 本 身 阶 次 高 结构 复 
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杂 , 使 其 直接 应 用 于 工程 实际 受到 很 大 的 阻碍 1. 
此 , 寻求 结构 简单 和 阶 次 低 且 具有 重 棒 控 制 性 能 的 
控制 器 , 已 成 为 鲁 棒 控制 发 展 的 一 个 新 的 研究 方 
向 9 

PID 控制 器 结构 简单 且 易 于 实现 , 迎合 了 人 们 
对 这 种 控制 器 的 强烈 需求 , 因而 在 工程 实践 中 得 到 
广泛 的 应 用 . 许多 学 者 对 工业 控制 过 程 中 这 类 参数 
不 确定 对 象 的 PID 控制 进行 了 长 期 的 研究 和 探讨 . 
Ziegler 和 Nichols 提出 了 著名 的 ZN 整定 公式 "1. 为 
克服 时 滞 现 象 对 控制 系统 的 影响 , 文献 [ 8] 提出 了 


时 滞 系 统 、PID 控制 的 研究 ; 盛 安 冬 ( 1964 一 ), 男 , 浙江 海 


盐 人 ;教授 ,博士 生 导 策 , 从 事 火 力 控制 ,满意 滤波 等 研究 . 
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著名 的 Smith 预 估 器 . 近年 来 , Silva 对 一 阶 参数 不 确 
定时 滞 系 统 的 鲁 棒 PID 控制 进行 深入 研究 , 推导 出 
求 取 所 有 和 鲁 棒 PID 控制 器 参数 区 域 的 算 湛 “4. 

一 般 情 况 下 , 工业 控制 过 程 中 的 系统 模型 参数 
存在 不 确定 性 , 而 以 时 澡 常 数 变化 最 为 显著 . 本 文 针 
对 时 滞 系 统 参数 不 确定 的 特点 , 给 出 了 和 鲁 棒 PID 控 
制 器 的 设计 方法 . 与 上 述 方法 相 比 , 本 文 算法 简单 ， 
并 可 迅速 得 出 所 有 重 棱 PID 控制 器 参数 集 , 从 而 提 
高 了 工业 控制 过 程 中 PID 控制 器 的 设计 效率 . 
2 问题 描述 

工业 控制 过 程 中 绝 大 多 数 系统 可 用 一 阶 时 滞 
系统 来 近似 ("1, 其 控制 系统 框图 如 图 1 所 示 . 其 中 : 
r( 为 系统 输入 ,u(t1) 为 控制 信号 ,y(t) 为 系统 输 
出 ; G(s) 为 参数 不 确定 一 阶 时 浪 系 统 , 即 


ee a Ee 


1 ”控制 系统 框图 


k Co 
Ts+ ie “， (J 


式 中 : E [有 1] 为 系统 静态 增益 ,2 ES [Li,L2| 为 
系统 时 滞 常 数 ,7 E [71,7,] 为 系统 时 间 常 数 ; C( s) 
为 控制 器 , 本 文 考 虑 PI 控制 器 和 PID 控制 器 两 种 情 


G(s) = 


形 , 即 
PI 控制 器 
C(s) = 而 + 1/s， (2) 
PID 控制 器 
C(s) = 而 + ki/s+ kas. (3) 


其 中 :为 比例 系数 ,i 为 积分 系数 ,fs 为 微分 系数 . 
如 何 整定 PI 或 PID 控制 器 参数 , 使 整个 对 象 族 G(s) 
稳定 , 是 本 文 所 要 解决 的 问题 . 
3 ”稳定 区 域 算法 
3. 1 PI 控制 情形 

对 于 式 (1) 描述 的 一 阶 时 湿 系 统 G(s) ,在 PI 控 
制 情形 下 , 闭环 控制 系统 的 特征 多 项 式 为 

6(s) = (Ts+ 1)s+ ke (ks+ i). (4 

使 s = jwWwE(0,+ co), 代 和 人 式 (4 ,得 

Oj = (jdor + Djw+ ke (knwW+ ki). (5) 
注意 到 e ”= cos(UL) - jsin( (7), 则 式 (5) 可 化 为 

0(jWw) = 

— WT + 有 nsin( OL) + kkicos(WL) + 

j[wW+ kkrWeos( WL) - kkisin( ul )]. (0) 
使 00 好 三 102 可 得 
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od + jicos(UL) = WT/E, 


Hpos(UD) -= pisin(UL) = - Wh. (7) 
解 方程 组 (7) , 得 
如 = [TOIn(W) - cos(OD)]/， i 


k= [TwWcos(WL) + Win(WL)]/E. 

根据 式 (8), 可 在 如 不 ;平面 上 绘 出 稳定 边界 线 
bp, i, 09 . 稳定 边界 线 将 心态 平面 分 割 成 稳定 区 
域 和 不 稳定 区 域 . 在 每 个 区 域 分 别 选取 一 个 测试 点 
(kb, i) , 则 可 确定 该 区 域 是 否 为 稳定 区 域 .绘制 稳定 
边界 曲线 时 , 频率 w 的 逼 历 取 值 非常 重要 . 根据 文献 
[9] 的 结果 , 今 *= jw wE(0,+ co), 当 Im C(Ojd 
= 0 时 , w 轴 被 方程 Im 6C(jd = 0 的 正 实 根 分 割 成 不 
同 的 区 间 , 测试 每 一 区 间 便 可 确定 是 否 存在 相应 的 
稳定 区 域 . 

记 w 为 方程 Im 6C(jd = 0 的 根 , 则 ww 必 满 足 
TW cos(W LIL) + sin( WL) = 0. 经 测试 可 知 , 若 7 
> 0, 则 只 有 当 w E(wm(2L), WL) 时 ,在 (0,wW ) 
上 存在 稳定 区 域 ; 若 7 < 0, 则 只 有 当 w € (0, 
(27)) 时 ,在 (0,w ) 上 存在 稳定 区 域 . 由 式 (8) 可 
得 请 (WU ,7Z) = 0, 则 w 为 如 ;平面 上 稳定 边界 线 第 
1 次 与 轴 相 交 时 所 对 应 的 w 值 (wj 友 0) . 解 方程 
Im G(jw) = 0, 可 求 得 , 即 只 要 在 ( 0,w ) 区 间 绘 
制 稳定 边界 线 , 便 可 得 出 所 有 控制 器 的 稳定 区 域 . 

当 被 控 系 统 G(s) 为 开 环 不 稳定 时 沿 过 程 时 , 使 
该 系统 稳定 的 PI 控制 器 的 存在 是 有 条 件 的 . 这 里 给 
出 PI 控制 器 的 存在 性 定理 . 为 表述 方便 ,将 式 (1) 改 
写 为 


he 
To le 
其 中 :LL > 0,k> 0,7 > 0. 
定理 1 对 于 式 (9) 描述 的 一 阶 不 稳定 时 滞 系 
统 , 存在 PI 控制 器 使 该 系统 稳定 的 必要 条 件 是 L/7 
<' |; 


证 明 根据 式 (8), 不 同 的 L/7 值 对 应 的 系统 在 


G(s) = (9) 


图 21| 不 同 L/ 研 伪 对 应 的 稳定 边界 线 
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名 帮 平面 上 的 稳定 区 域 如 图 2 所 示 . 由 图 2 可 知 , 若 
稳定 的 PI 控制 器 存在 , 则 有 


人 (10) 

即 

jim OE - 

oo dho 

jm EY) TEsin( WF) - sin( WF) - Wcos( WL) _ 
a Tsin(wDJ + TLavos(oL) — Lsin( CL) 

27 sin( ( 忆 sin( WL) 

a T SO TL i CD — cos(UuL) 

w>0+ T 于 Toos(W)- L ue 

2(1- 了 /了 

Ro (11) 


由 此 可 得 L/T < 1 或 L/7 > 2. 由 图 2 知 L/7 > 2， 
不 存在 稳定 区 域 , 故 舍 去 , 从 而 得 L/T < 1. 
3. 2 PID 控制 情形 
与 PI 控制 情形 类 似 , 对 于 式 ( 1) 描述 的 一 阶 时 
滞 系 统 , 对 应 式 (7) 不 难得 到 
jin((L) + Hicos(() — 
kawcos( UL ) = WT/E, 
kWeos (WL) -= kisn(WOL) + 
hawsin( OL) = — wk. 
固定 值 , 可 解 得 
如 = [TOIn(WL) - cos(0W)1/Eh, 
k= [Tujcos(0L) + usin(WL) + (13) 
kkatd | /Ek. 
根据 式 (13) ,在 hoi 平面 上 绘 出 稳定 边界 线 ， 
便 可 得 出 对 应 的 稳定 区 域 . 注意 到 式 ( 13) 中 心 值 与 
ja 值 无 关 , 故 可 通过 瑶 历 hi 值 , 在 如 大: 平面 上 绘 
稳定 边界 线 , 以 确定 所 有 稳定 的 值 区 间 . 同 理 ， 
定语 值 , 在 刀 趟 (平面 上 亦 可 确定 对 应 的 PID 控制 器 
稳定 区 域 
当 固 定 思 值 时 , 解 方程 (12) 存在 分 尽 为 雾 的 情 
形 , 故 不 能 由 此 方法 确定 kki 平面 上 PID 控制 器 的 
稳定 区 域 . 对 同一 心 值 ,大 平面 上 PIDD 控制 器 的 稳 
定 区 域 是 由 若干 条 坦 线 围 成 的 凸 集 t". 根据 这 一 特 
点 ,结合 馈 五 ; 平面 和 4 平面 上 的 稳定 边界 线 , 选 
取 对 应 的 交点 , 同样 可 确定 ha; 平面 上 PID 控制 器 
的 稳定 区 域 . 
当 被 控 系 统 G(s) 为 开 环 不 稳定 时 滞 过 程 时 ， 
PID 控制 器 的 存在 是 有 条 件 的 . 这 里 给 出 PID 控制 
器 的 存在 性 定理 : 
定理 2 ”对 于 式 (9) 描述 的 一 阶 不 稳定 时 滞 系 
统 , 存 在 Pip 控制 器 使 该 腔 统 稳定 的 必要 荣 件 是 
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天 了 和 2 
注意 到 式 (13) 的 hs 值 与 如 值 无 关 , 故 由 ow 
平面 上 的 稳定 边界 线 , 仿照 定理 1 即 可 证 明 本 定理 . 
4 和 鲁 棒 PID 控制 器 设计 
n 阶 时 澡 系 统 的 叭 特征 多 项 式 可 描述 为 
p(s)= aoos 二 2 [名 Ce Ws "(14 


其 中 :qi = Qt+ jBi, Oi, Bi ER(i Gz 0) 为 常量 , ao 
天 0:0= ho < hi< …< hy,h 为 对 应 于 系统 时 洁 


we 


浓 数 的 常量 . 

对 于 n 阶 参数 不 确定 时 语系 统 , 其 准 特 征 多 项 
式 族 可 描述 为 

PpP*. 


n N 
(p(s) p(s) = ao0s + > ( 2 ire 3 


(aoo a10, ,QIN ,anN) EF,aw GA 0)， ( 19) 
其 中 EC” ”为 不 确定 参数 集 . 本 文 考 虑 的 参数 
不 确定 时 灌 系 统 族 准 多 项 式 族 , 是 由 n 阶 准 多 项 式 
Pi1(s) ,Pp2(s),…,pr(s) 所 生成 的 准 多 项 式 多 面体 , 即 

P *» conv{p1(s),p2(s),**,p 7(s)}, (10) 
其 中 conv 表示 山 生 成 .着 P 为 实 多 项 式 族 , 则 准 多 
项 式 p(s) 首 系数 必须 同 号 

定义 并 ” ” 若 D 为 复 平面 上 对 称 于 实数 轴 的 
任 一 区 域 , 则 称 时 滞 系 统 (1) 是 D 稳定 的 , 当 且 仅 当 
由 式 ( 14) 描述 的 叭 多项式 的 所 有 零点 均 在 九 内 . 若 
D 为 左 半 平 面 , 且 由 式 (14) 描述 的 准 多 项 式 的 所 有 
零点 均 在 D 内 , 则 称 p (s) 是 稳定 的 . 

引 理 1 考虑 由 式 (16) 描述 的 n 阶 准 多 项 式 
多 面体 P,D 为 复 平面 上 对 称 于 实数 轴 的 任 一 区 域 ， 
x 和 y 为 复 平面 上 的 两 点 , 满足 对 任 一 给 定点 x E 
D"(D" 表示 DD 在 复 平面 的 补 集 ) 和 任 一 M > 0, 存在 
实数 Q. 若 点 y 满足 y 二 MM 且 Re(y) 二 %% 则 在 D"“ 
内 总 能 找到 一 条 由 x 通 向 y 的 连续 路 径 , 使 P 是 D 
稳定 的 , 当 且 仪 当 P 的 所 有 楼 边 均 是 D 稳定 的 . 

对 于 式 (1) 描述 的 一 阶 时 洁 系 统 , 倒 


hk; a 
7 Ls . 。 a 
Gi(s) 二 Ts 1° 1, 二 1, 分 站 三 1 , 2， ,8, 


(17) 
则 一 阶 参数 不 确定 时 滞 系 统 族 可 描述 为 
T(s) = 2 NGi(s), NX 二 0, 和 1,r= 8. 
(18) 
定理 3 ”考虑 式 (18) 描述 的 一 阶 参数 不 确定 


时 灌 采 统 族 了 ( 洲 用 PI 入 制 器 进行 控制 ,使 恋 采 统 
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族 7(s) 稳定 , 当 且 仅 当 用 同一 PI 控制 器 使 每 一 


Gi(s) 稳定 . 
证 明 “” 兮 闭环 控制 系统 的 开 环 传递 本 数 为 
(0 EC , rb 
Fi(s) Di) 1,2,…%,r= 8, (19) 


其 中 pi(s) 和 gi(s) 为 准 多 项 式 . 则 由 式 ( 18) 描述 的 
一 阶 参数 不 确定 时 涨 系统 族 的 特征 多 项 式 可 描述 为 

P * conv{(Fi(s) + 1)p1(s), (F2(s) + 

Dp2s), ,Fr(s) + Dp"(s)}. (20) 

取 D 为 复 平面 左 半 平 面 ,对 任 一 给 定点 x ED“ 和 任 
一 M > 0, 存在 实数 若 点 y 满足 vy 二 M 且 
Re(y) 之 % 则 在 Dp" 内 总 能 找到 一 条 由 x 通 向 y 的 
连续 路 径 . 根据 引 理 1, 要 使 7(;s) 稳定 , 只 需 使 7( 39) 
的 准 特 征 多 项 式 族 P 的 每 一 楼 边 (Fi(s) + 1)pi(3) 
是 D 稳定 的 即 可 . 注意 到 (Fi(s) + 1)pi(s) 为 Gi(3) 
的 闭环 特征 多 项 式 , 使 每 一 ( Fi( s) + 1)pi(s) 是 D 稳 
定 的 即使 Gi(s) 是 稳定 的 . 

定理 4 考虑 式 (18) 描述 的 一 阶 参数 不 确定 
时 洁 系 统 族 7T(s), 用 PID 控制 器 进行 控制 , 使 该 系 
统 族 7(s) 稳定 , 当 且 仪 当 用 同一 PID 控制 器 使 每 一 
Gi(s) 稳定 . 

仿照 定理 3 的 证 明 很 容易 证 得 本 定理 . 

对 于 实际 工程 控制 器 设计 而 车 ,满意 控制 ” 认 
为 :控制 系统 设计 要 求 不 是 某 单一 性 能 指标 达到 最 
优 , 而 是 必须 同时 满足 由 多 种 性 能 指标 所 构成 的 控 
制 性 能 期 望 指标 集 . 根据 实际 工程 对 控制 器 设计 的 
要 求 , 可 选择 相应 的 性 能 指标 , 如 动态 误差 系数 、 系 
统 最 大 开 环 增益 、 调节 时 间 、 超 调 量 等 时 域 性 能 指 
标 , 以 及 幅 值 众 度 、 相 角 容 度 等 开 环 频 域 性 能 指标 和 
和 鲁 棒 性 能 指标 , 在 所 求 PID 控制 器 参数 区 域内 选择 
同时 满足 各 种 性 能 指标 的 参数 , 便 可 使 被 控 系 统 达 
到 所 期 望 的 性 能 要 求 . 
5 数值 算 例 

例 1 考虑 文献 [13] 中 一 阶 参数 不 确定 时 灌 


Pog 
ke 


Ts+ 1’ 
其 中 :LE[1,2],k 全 [1,2],7 E[1,2]. 由 式 (17) 
可 得 如 下 8 个 楼 边 系统 : 


G(s) = (2 
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李 银 份 等 : 参数 不 确定 时 注 系 统 的 鲁 棒 PID 控制 
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根据 3. 1 节 PI 控 制 器 稳定 区 域 算 法 ,bi; 平 面 
上 8 个 棱 边 系统 所 对 应 的 稳定 区 域 如 图 3 所 示 . 由 
定理 3 可 知 , 当 且 仅 当 8 个 棱 边 系统 同时 稳定 时 , 则 
G(s) 即 为 稳定 的 . 图 3 中 阴影 部 分 为 8 个 棱 边 系统 
稳定 区 域 的 交集 , 即 为 所 求 PI 控制 器 的 稳定 区 域 . 
在 该 区 域 上 任 选 控制 器 参数 (种 , 如 ) , 系统 (21) 均 为 
稳定 的 . 
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图 3 ” 例 1 PI 控制 器 的 稳定 区 域 
例 2 考虑 如 下 一 阶 参 数 不 确 定时 滞 过 程 : 
ke” 

Ts— 1 
其 中 :LE [1,1.2],k [0.8,1],7T Ef[l,1.5]. 则 
存在 L/7T > 1 的 情况 , 由 定理 1 知 ,不 存在 PI 控制 器 
使 该 系统 族 稳定 ;但 L/T < 2, 故 可 用 PID 控制 器 进 
行 控制 . 由 式 (17) 可 得 如 下 8 个 楼 边 系统 : 


G(s) = (22) 


Gt) eT Gals) = ee 1 
人 运 Ps) = 工 5 1 
0. 8e 1* 0.8e '” 
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图 4 例 2 PID 控制 器 的 稳定 区 域 
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根据 3.2 节 PID 控制 器 稳定 区 域 算法 , 取 f4 = 
0. 75, 局 三 平面 上 8 个 楼 边 系统 所 对 应 的 稳定 区 域 
如 图 4 所 示 . 由 定理 4 知 , 8 个 棱 边 系统 稳定 区 域 的 
集 即 为 所 求 控 制 器 的 参数 区 域 . 在 该 区 域 上 任 选 


控制 器 参数 (, hi, ka = 0.75), 系统 (22) 均 为 稳定 
的 . 
6 结 论 


本 文 给 出 一 种 快速 计算 所 有 和 鲁 棒 PI 和 和 鲁 棒 
PID 控制 器 参数 集 的 方法 . 该 方法 首先 基于 推广 到 
时 浊 系 统 的 楼 边 定理 , 在 kr; 平面 上 绘制 各 个 楼 边 
系统 的 稳定 边界 线 ; 然后 确定 所 有 稳定 区 域 的 交集 ， 
即 为 所 求 控制 器 参数 集 . 数值 算 例 表明 了 该 方法 的 
优越 性 和 有 效 性 . 
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